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I. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length 

l 

t 

F 

second - . - 

m 

s 

kg 

foot (or mile) 

second (or hour) _ = — 

ft (or mi) 
sec (or hr) 
lb 

Force. 

weight of 1 kilogram 

weight of 1: pound 

p„„,. 

P 

V 



horsepower 

fps 

Speed 

/kilometers per hour.. 

(.meters per second 

kpli 

mps 

miles 'per hour 

feet per second 


2. GENERAL SYMBOLS 

Weight = mgr v Kinematic viscosity 

Standard acceleration of gravity=9.80665 m/s 2 p Density (mass per unit volume) 

or 32.1740 ft/sec 2 Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at. 15° C 

W and 760 mm; or 0.002378 lb-ft -4 sec 2 

Mass=-— Specific weight of "standard” air, 1.2255 kg/m 2 or 

Moment of inertia =m/r. (Indicate axis of 0.07651 lb/cu ft 
radius of gyration k by proper subscript.) 

Coefficient of viscosity 

3. AERODYNAMIC SYMBOLS 


Area 

.Area of wing 
Gap 
Span 
Chord 


Aspect ratio, g- 
True air speed 
Dynamic pressure, -^pV* 

Lift, absolute coefficient Cx=-4s 


([S 

D 


Drag, absolute coefficient 6 r D =^ 

Profile drag, absolute coefficient C Dq ~’^ 


D n 


Di 


Induced drag, absolute coefficient (7 D 
Parasite drag, absolute coefficient C' D „='^| 


C 


Cross-wind force,- absolute coefficient Cc-=^, 


Q 

Q 

R 


a 

€ 

«< 

a. 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 
VI 

Reynolds number, p— where l is a linear dimen- 
sion (e.g., for an airfoil of 1 .0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured fiom zero- 
lift position) 

Flight-path angle 
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THEORETICAL SYMMETRIC SPAN LOADING AT SUBSONIC SPEEDS 
FOR WINGS HAVING ARBITRARY PLAN FORM 

By John DeYoung and Charles W. Harper 


SUMMARY 

A method is shown by which the symmetric span loading for 
a certain class of wings can be simply found. The geometry 
of these wings is limited only to the extent that they must have 
symmetry about the root chord, must have a straight quarter- 
chord line over the semispan, and must have no discontinuities 
in twist. A procedure is shown for finding the lift-curve slope, 
pitching moment, center of lift, and induced drag from the 
span load distribution. A method of accounting for the effects 
of Mach number and for changes in section lift-curve slope is 
also given. 

Charts are presented which give directly the characteristics 
of many wings. Other charts are presented which reduce the 
problem of finding the symmetric loading on all wings falling 
within the prescribed limits to the solution of not more than four 
simultaneous equations. 

The loadings and wing characteristics predicted by the 
theory are compared to those given by other theories and by 
experiment. It is concluded that the results given by the subject 
theory are satisfactory . 

The theory is applied to a number of wings to exhibit the 
effects of such variables as sweep, aspect ratio, taper, and 
twist. The results are compared and conclusions drawn as to 
the relative effects of these variables. 

INTRODUCTION 

Theory predicts and experiment has shown that the effects 
of compressibility on wing characteristics can be delayed and 
diminished through the use of wing sweep and/or reduced 
aspect ratio. Experiment has also shown that if sufficient 
wing sweep or reduction in aspect ratio is used to redpce 
compressibility effects, then the characteristics of the wings 
will be different from those of conventional wings. The 
importance of being able to predict the characteristics of 
these heretofore unstudied plan forms in incompressible 
flow is obvious. Not only does it become possible to make 
a systematic study of the effects of major plan-form changes, 
but also a base is supplied from which the subcritical effects 
of compressibility can be predicted. 

Lifting-line theory in the past has been so modified and 
extended that the characteristics of wings having no sweep, 
moderate to high aspect ratio, and any taper ratio can be 
determined readily with good accuracy. For studies re- 
juiring a higher degree of accuracy, lifting-surface theories 


have been used, but generally it has been found that the 
additional complexity of these methods has not sufficiently 
improved the predictions to warrant common use. 

Lifting-line theory, however, proved wholly inadequate 
when used to predict the characteristics of wings having 
appreciable angles of sweep and/or very low aspect ratio. 
Lifting-surface theories, in contrast, made satisfactory pre- 
dictions of the characteristics of these wings although the 
extent of the computing labor involved prevented the under- 
taking of any general study. 

Continuing studies of the problem resulted in the theory 
first satisfactorily presented by Weissinger (reference 1) for 
the case of wings having a straight quarter-chord line across 
the wing semispan. It was found that this approach, which 
can be considered a simplified lifting-surface theory, enables 
rapid and satisfactory predictions to be made of the incom- 
pressible-flow characteristics of wings having swept and/or 
low-aspect-ratio plan forms as well as those of more con- 
ventional plan form. Further, it was found to be admirably 
suited for the problem of undertaking a systematic study of 
the effects of plan form on wing characteristics. Develop- 
ment of the method and procedures for its application to- 
gether with the results of applying it to study the charac- 
teristics of a series of wings have been presented in references 
2, 3, and 4. It is the purpose of this report to combine and 
extend the material contained in these three references. 

NOTATION 

m number of span stations at which circulation and 

downwash are found 

n an integer defining a spanwise station on the wing 

quarter-chord line for which the value of circula- 
tion is determined 

v an integer defining a specific point within the wing 

plan form for which the boundary condition of 
no flow through the wing is applied 
V dimensionless lateral coordinate measured per- 

pendicular to the plane of symmetry, fraction of 
semispan 

Vc.p. spanwise location of center of pressure 

V c.p- a spanwise location of center of pressure due to addi- 

tional lift on the wing 

<t> trigonometric spanwise coordinate (cos- 1 f), radians 
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a wing angle of attack measured in a plane parallel 

to the plane of symmetry, degrees or radians 
a, angle Of attack at station v, degrees or radians 

«' 0 angle of attack at station v for zero net lift on the 

wing, degrees or radians 

a , 0 angle of attack of the wing-root section for zero net 

lift on the wing, degrees or radians 
e, angle of twist of the wing at station v measured 

relative to the wing-root section («»=a» 0 — <*r 0 ), 
degrees or radians 

G n dimensionless circulation identical to the load 
coefficient at spanwise station n 


G„ 0 value of G„ for zero net lift on the wing 

Gn a value of G„ due to additional lift on the wing 

A y „ coefficient dependent on wing geometry and in- 
dicating the influence of arbitrary loading at span 
station n on the downwash angle at span station v 
a yn coefficient dependent on wing geometry and in- 

dicating the influence of symmetric loading at 
span station n on the downwash angle at span 
station v 

b span of the wing measured perpendicular to the 

plane of symmetry, feet 

c local chord of the wing measured parallel to the 

plane of symmetry, feet 
c, local chord through point v, feet 


A 

S 

A 


Ci 

c i a 

c ‘o 

Cl„ 


C L 

Ol. 

Cm 

c mo 

C m . . 


/J>\ 

\L cdv ) 


feet 


c a » average chord f -g y feet 

M.A.C. mean aerodynamic chord 

X taper ratio 

aspect ratio ^ s) 
wing area, square feet 

geometric sweep angle of the wing quarter-chord 
line, positive when quarter-chord line is swept 
aft of a line normal to the plane of symmetry, 
degrees 

Ag compressible sweep parameter (^Ag= tan -1 - ^ 


section lift coefficient 


(I) 


section lift-curve slope, per radian or per degree 
section lift coefficient for zero net lift on the wing 
section lift coefficient due to additional lift on the 
wing 

wing lift coefficient 


(§) 


lift-curve slope, per radian or per degree 

pitching-moment coefficient xO 

pitching-moment coefficient for zero net wing lift 
on the wing 

pitching-moment coefficient due to additional lift 
on the wing 



a.c. 


r 

h 

p 

w 


V 

<L 

M 


Ky 

Hy 

dy 


induced-drag coefficient ^ i n d u g pd drag ^ 

induced-drag coefficient due to basic loading (zero 
net lift on the wing) 

induced-drag coefficient due to additional loading 
(net lift on the wing) 

longitudinal position of aerodynamic center, meas- 
ured from the leading edge of the mean aero- 
dynamic chord, in percent of the mean aero- 
dynamic chord 

circulation, square feet per second 
absolute distance from a vortex to a downwash 
point, measured perpendicular to the vortex, feet 
density of air, slugs per cubic foot 
induced velocity, perpendicular to the mean chord 
line of the wing, positive for downwash, feet per 
second 

free-stream velocity, feet per second 
dynamic pressure, pounds per square foot 
compressibility parameter (-\ /l—M 2 ) 

Mach number 

ratio of the experimental section lift-curve slope Ci a 

2ir 

to the theoretical value of both taken at the 

same Mach number 
value of k at spanwise station v 


wing geometric parameter d 



scale factor 


METHOD 

THEORY AND LIMITATIONS 

A detailed mathematical development of the method is 
given in the appendix to this report. However, as a pre- 
liminary to the discussion of the use of the method and the 
results obtained from it, an outline of the theory and its 
application is given here. 

The wing is replaced by a plate of zero thickness but having 
a plan form and twist identical to that of the wing. It is 
assumed that the chordwise distribution of loading on the 
plate can be concentrated into a lifting line lying along the 
quarter-chord line of the wing. The method as developed 
here requires that the lifting line and, hence, the quarter- 
chord line of the wing in question, be a straight line over 
each semispan. At the plane of symmetry the method 
allows an angular discontinuity of the quarter-chord line to 
exist which therefore enables consideration of any degree of 
wing sweep . 1 The boundary condition which fixes the span- 
wise variation of the circulation is that the slope induced in 
the flow field by the downwash normal to the plate due to 
the lifting line and its trailing vortex sheet shall be equal to 
the slope of the plate with respect to the free stream at 
specified points (control points) within the wing boundary, 
or, in effect, that no flow shall occur through the plate at 
the control points. On the basis of two-dimensional theory, 
the chordwise location of these control points is chosen to 


1 It should be noted that this specifies that the angle of sweep is taken with respect to the 
quarter-chord line. In this report it will be understood that sweep angle refers to sweep of 
the quarter-chord line. 
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be at the three-quarter-chord line. 2 As will be shown later, 
this implicitly specifies that for each section the lift coefficient 
increases at the theoretical two-dimensional rate of 2r per 
radian change in true angle of attack of the section. The 
spanwise locations of the control points are chosen for 
reasons of mathematical expediency alone since only at 
certain fractions of the wing semispan is it convenient to 
solve the integral expressions for the total downwash. 

The method treats the lifting line and its trailing vortex 
sheet as being continuous. However, explicit values for the 
circulation strength are found only for those spanwise points 
on the quarter-chord line which correspond to the spanwise 
location of the control points. The total number of control 
points m can be as great as desired. Comparison of the 
loadings predicted where various values of m are used shows 
that little change occurs where m is made greater than seven. 
Further, comparison with experimental results has shown 
that as few as seven (one at midspan and six distributed 
symmetrically about this point) will enable good predictions 
to be made of the span loading of wings without sharp dis- 
continuities in span loading as would result from partial- 
span flaps. All of the material presented in this report is 
based on the seven-point solution with the exception of the 
appendix which is not restricted. 

Using the vortex pattern substituted for the wing and 
applying the boundary condition allows a set of simultane- 
ous equations to be formed, each of which involves (1) the 

slope of the flow ^ (or angle of attack of the plate 3 a,) 
induced by the total downwash at each control point v, (2) 
the load coefficient of the lifting line at each span- 

wise point n on the quarter-chord line, and (3) influence 
coefficients A, n which relate the influence of the circulation 
at any point n to the downwash at any point v and are a 
function of wing geometry only. The method shows that 
for an arbitrary loading the equations have the following 
form: 

( w \ J 71 . 

y ) =ot t =22Ar n Q n , r= 1,2, . . . m (1) 

* It is important for the reader to realize here that a choice has been made between a number 
of possible procedures. These possibilities arise from the fact that the exact location, on a 
tapered wing, of constant-percent- chord lines depends upon the orientation of the reference 
line along which the chord is measured. The orientation of the reference line is usually 
chosen such that an airfoil section so defined will have aerodynamic characteristics closely 
resembling those found two-dimensionally for the same section. This then enables an esti- 
mation of the effect of changes in section characteristics on over -all wing characteristics. For 
unswept wings, there is little reason to consider any orientation of the reference line other 
than parallel to the free-stream direction. However, when a wing is swept, the question of 
orientation of the reference chord cannot be so easily answered. Insufficient experimental 
data exists to determine the most satisfactory orientation, and strong arguments can be 
presented for at least two orientations, namely, parallel to the free stream .and perpendicular 
to some swept reference line. In the present analysis the reference chord was chosen as being 
parallel to the free stream since it greatly simplifies the mathematical procedure and since 
consideration of the differences expected from use of the alternate choice indicates they will 
be small, 

* The reader should note that the boundary condition is given by w ¥ = V , sin from 
which (*). is seen equal to sin a.. The substitution of a , for sin a, has the effect of 
increasing the value of loading on the wing above that necessary to satisfy the boundary 
condition. However, the boundary condition was fixed assuming that the shed vortices 
moved downstream in the extended chord plane. A more realistic picture is obtained if the 
vortices are assumed to move downstream in a horizontal plane from the wing trailing edge. 
It can be readily seen that, if this occurs, the normal component of velocity induced by the 
trails at the three-quarter-chord line is reduced and, if the boundary condition is to continue 
to be satisfied, the strength of the bound vortex must increase. It follows that substitution 
of av for sin a» then has the effect of accounting for the bend up of the trailing vortices. It is 
not known how exact the correction is, but calculations and experimental verification show it 
of to be the correct order. 


3 

Each equation gives the downwash angle at the control 
point v, the spanwise location of which is defined by 

vi r 

»? = COS -g- 

where the downwash results from the circulation at m points 
n on the wing the spanwise locations of which are defined 
also by 

nr 

77= cos -g- 

In the case of symmetrically loaded wings, each panel pro- 
duces an identical equation for the corresponding semispan 
point. Since only one of these identical equations is of 
value, the total series reduces to the equations correspond- 
ing to the wing midpoint and one panel,. For the seven- 
point solution, equation (1) is therefore written 

&r— y j (i„ n G n , v=l,2,3,4 (2) 

72 = 1 

where a vri represents the influence coefficients for the sym- 
metric seven-point solution. The set of four simultaneous 
equations so formed can be easily solved to obtain the dis- 
tribution of total load (in terms of G n ) on any wing for which 
the angle of attack at each spanwise station, sweep, and 
chord distribution are specified. The distribution of load 
is specified at only four spanwise stations, namely tj= 0.924, 
0.707, 0.383, and 0 {n— 1, 2, 3, and 4, respectively). Values 
of loading at additional spanwise stations can be found by 
means of the interpolation function given in the appendix 
(equation (A52)). 

The simplicity of the procedure depends to a large extent 
on the fact that the solution can be found in terms of the 
coefficients a, n . Even where these must be computed for 
each wing plan form the method offers computational ad- 
vantages over other equally accurate methods. However, 
because these a rn coefficients are a function of geometry alone, 
it is possible to relate them in a simple manner such that a 
limited amount of computation will give the a yn coefficients 
for all plan forms to which the method is applicable. Details 
of this procedure and the results of applying it are discussed 
in a later section of the report. 

The method assumes that the flow follows the wing surface 
and makes some allowance for the trailing sheet aft of the 
trailing edge becoming horizontal. 3 Hence, the method 
should apply to higher angles of attack with considerable 
accuracy, provided the flow remains along the wing surface. 
The method assumes incompressible flow but it will be shown 
how the effects of compressibility can be included within the 
limits of applicability of the Prandtl-Glauert rule. The 
method assumes the theoretical section lift-curve slope 0/ 
2ir (or with account taken of compressibility, 2x//3) but a pro- 
cedure will be shown which accounts for the variation in 
section lift-curve slope from the theoretical value. 

It is clear from the foregoing outline of the theory that the 
method can account for varia tions in those geometrical char- 

* See footnote 5. 
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acteristics of wings, namely plan form and twist, which have 
the greatest influence on the spanwise distribution of lift. 
With the exception of variations in c !a the method cannot 
directly account for any effects dependent upon the geometry 
of the airfoil section even though these may affect the span 
loading. The substitution of uniformly cambered sections 
for uncambered sections across a wing span is assumed to 
change only the wing angle for zero lift and this change is 
assumed equal to that shown by the section. The substi- 
tution of variable camber is assumed equivalent to twisting 
uncambered sections. 

PROCEDURE FOR DETERMINATION OF AERODYNAMIC CHARACTERISTICS 
FROM SPAN LOADING 

The foregoing section has shown a method by which the 
symmetric span load distribution of any wing having a 
straight quarter-chord line over the semispan can be deter- 
mined from a knowledge of wing geometry only. With 
such loading determined it becomes possible to quickly find 
other characteristics. It should be noted, however, that 
these characteristics are derived directly from the wing load 
distribution and that no further aerodynamic theory is in- 
volved. 

It is possible to find the gross load distribution and re- 
sultant characteristics directly for a wing at any angle of 
attack and having any plan form and twist. Past experience, 
however, has shown that gross characteristics can better be 
studied if the basic and additional type loadings are handled 
separately. Since the two types of loading are additive, 
this procedure is permissible. 

Basic loading is that existing with zero net lift oh the wing 
and is, therefore, due to twist or effective twist 4 (e. g., span- 
wise change in camber) of the wing-chord plane. The basic 
loading and characteristics dependent on it are unchanged by 
the addition of load due to uniform spanwise wing angle- 
of-attack change and are equal at all angles of attack to that 
found for zero net lift on the wing. 

Additional loading is that due to equal geometric angle-of- 
attack change at each section of the wing. The distribution 
of additional load is a function only of wing plan form and is 
thus independent of any basic load due to twist existing on 
the wing. The magnitude of the additional load is a function 
only of angle of attack of the wing and thus each equal 
increment of angle of attack will give the same increase and 
distribution of additional load irrespective of the gross load 
on the wing. The wing characteristics due to additional 
load of any given plan form are thus a function of the lift 
coefficient or angle of attack of the wing. 

In the following sections the procedure for determining 
basic-type span load distribution and the characteristics 
associated with it (denoted by a subscript 0) is first presented 
and then the procedure for finding additional-type span load 
distribution and the associated characteristics (denoted by a 
subscript a). Finally, it is shown how characteristics due 
to gross load distribution (denoted by absence of a subscript) 
can be found. 


4 Hereafter reference will be made to twist only. The reader will understand that effective 
twist will be bandied in an identical manner. 


PROCEDURE FOR DETERMINATION OF AERODYNAMIC CHARACTERISTICS 
DUE TO BASIC LOADING 

Span load distribution and angle of zero lift for arbitrary 
twist. — It is not possible to obtain the distribution of basic 
load on a twisted wing directly from equation (2) since the 
values of a, for each station are not generally known for the 
condition of zero net lift and thus, as written, eight unknown 
values appear. However, only one unknown has actually 
been added to the four unknowns (the individual loads of 
equation (2)) since, while four values of a„ appear, they are 
related to one another through knowlege of the twist distri- 
bution. To form the fifth equation required for the solution 
in addition to the four represented by equation (2), use is 
made of equation (A46) of the appendix which gives the 
total lift in terms of the four individual loads. Thus 

(<? 4 + 2 S G n sin <*>„) (3) 

Equating this to zero (basic loading) and with equation (2), 
a set of five simultaneous equations is formed, the solution of 
which will give the values of loading at the spanwise points 
and the angle of zero lift of the reference station. Thus 

3 

0= 6*4+2 s Gn sin 4> n 

n~\ 

4 

I ^refn ^ ; &vnG „, V 1 , 2, 3, 4 

71 = 1 

Solution of the five equations will give values of the angle of 
zero lift of the reference section and the four individual loads. 
The angle of attack of sections other than the reference sec- 
tion is determined from the twist. 

Some simplification of this process can be made, however. 
When equated to zero lift, equation (A46) involves only the 
four individual loads, and it is therefore possible to express 
one load in terms of the other three. Elimination of one 
unknown in the remaining equations enables a solution to 
be made for the angle of zero lift and three individual loads 
from a set of four simultaneous equations. Finally, since 
the individual section angles for zero lift are each expressed 
as the sum of the reference angle and a twist angle (which is 
zero for the reference station), one equation and one un- 
known (the reference angle) can be eliminated by subtracting 
one equation from the remaining three. Values can then be 
found for the three loads from solution of only three simul- 
taneous equations, with the angle of zero lift for the refer- 
ence chord and the value of the fourth load found from the 
previously eliminated equations. 

The exact form of the equations depends on the section 
chosen as a reference. If, as is customary, the root is taken 
as a reference, then 


y ( 4 ) 


a Vo =a 4o +e,=a ro +e, (5) 

where e, is the angle of twist at station v with respect to 
station 4, and the set of equations required for the solution 
take the following form. The three simultaneous equations 
giving the three outboard loads are 
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<1 = [«U — *41 — 0.765 (*14 — fl44)]^l 0 +[®12 — *42 — 1 .414 (*14 — *44)]f*2 0 +[*13 — *43 — 1.848 (*14 — *44)]f?3 0 
<2 = [*2l — *41 — 0.765 (*24 — ®44) ]f?l 0 + [*22 *42 1.414 (d 2 4 — <*44) ]f?2 # + [®23 — *43 — 1 • 848 (ffl 2 4 — *44) ] > 

e 3 = [(I 31 — a 41 — 0.765 ( 4 X 34 — * 44 ) ] [®32 *42 1.414 (<134 O 44 ) ] f? 2 0 + [*33 — *43 — 1-848 (034 — * 44 ) ] G 3q 
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( 6 ) 


The equation giving the root load is 

G t = - (0.765Gf Io +1.414Gf 2o +1.848<?3 0 ) (7) 

The equation giving the angle of the root section for zero 
lift is 

OV 0 =a4 0 = (4X41 — 0 . 7 65*44) $1 0 + (*42 — 1.414*44) Gz a 

+ (*43 — 1 . 848 * 44 ) <?3 0 (8) 

The dimensionless circulation G n can be related to other' 
forms of loading coefficients by 


of the appendix) and summing the moment of each wing 
element about an arbitrary axis, the value of the couple due 
to the basic load distribution can be found. About an axis 
normal to the plane of symmetry and in terms of dimen- 
sionless circulation this is 

Ah 

Cm <‘ == ~ MAG tan A (°- 138 Cf >o+ 0 - 198 ^o+°- 135 ^o+ 
0 . 016 <? 4 0 ) (10) 

If equation (4) is again used to eliminate G i(1 , equation (10) 
becomes 



Conversely, if the load at each of the four span stations is 
known for the case of zero wing lift, it becomes a simple 
matter to solve equation (2) directly and determine the 
required wing twist. 

Induced drag. — With the values of loading found at the 
four spanwise positions, the induced drag can be found from 
the following equation which is derived in the appendix. 
In terms of the nondimensional circulation G n , 

C Dio =^[G\+G\+G\+^-G io (OmQG lo + 

0.789G 3o ) — £2 0 (0.733<7 1 o + 0.845£ 3o )] (9) 

Through use of equation (4), the induced drag can be put 
in terms of only three values of G„. Thus 

C D<o =;rA[0.6686 !2 i 0 +G !2 2 0 +2.082G !2 3 0 -t- (0.2156^+ 

1 . 442(t 3o ) f? 2 0 + 1 .06 1 G lfj G 3o ] (9a) 

Pitching moment.- — It is quite evident that, as a wing is 
swept, the distribution of span load will increasingly affect 
the wing pitching moment where this is measured about an 
axis normal to the plane of symmetry. Correspondingly, as 
the wing is swept, any changes in the chord load distribution 
exert relatively less effect on wing pitching moments. Fur- 
ther, since the possible distortion of a chord loading from the 
conventional is not great, the span load distribution begins 
to control the wing pitching moments even at small angles 
of sweep. Since the subject method gives, as will be shown, 
good predictions of span load distribution, then it is safe to 
assume that good estimations of pitching moments should 
be possible. 

The moment produced by basic-type load distribution is 
in the form of a pure couple since zero net load exists on the 
wing. It follows immediately that the magnitude and sign 
of the moment are independent of the fore-and-aft location 
of any reference axis normal to the plane of symmetry. 
Using the expression for continuous loading (equation (A 17) 


Cm ^ = ~ MAG tan A (°- 126(? io+ 0 - 175 ^o+ 0 - 1066r 3 0 ) 

(10a) 

PROCEDURE FOR DETERMINATION OF AERODYNAMIC CHARACTERISTICS 
DUE TO ADDITIONAL LOADING 


Span load distribution and lift-curve slope. — Since the 
additional loading distribution for wings of a given plan 
form is the same, regardless of twist or camber, it is necessary 
to consider only the case of the flat wing. Further, since 
the magnitude of the additional load varies directly with 
angle of attack, it is only necessary to consider unit angle- 
of-attack changes as all other loadings will be directly pro- 
portional to this. Thus, to find the additional load distri- 
bution, equation (2) is written in the following form 

l = r=l,2,3,4 (11) 

71=1 a 


Solution of the set of four simultaneous equations gives the 

G n 

values of circulation per radian ~ at the four spanwise 

stations t;= 0.924, 0.707, 0.383, and 0. 

G n 

Substitutions of the values of — - in the expression for lift 

coefficient gives the wing lift coefficient for one radian 
change in angle of attack, or in effect, lift-curve slope. 
Thus, with equation (3) 


dCL =G = *A(G± 


da 


8 V 


-+1.848 — +1.414 —^+0.765 


( 12 ) 


Gn a 

The dimensionless circulation per radian — can be ex- 
pressed in the more usual loading-coefficient form through 
the relations 


G n«_ G L «( C 'a C \ _ r f C ‘a C \ 

a - 2 A \C L cJ n \2 bC L Jn 

Induced drag. — The induced drag due to additional load- 
ing can be computed at any lift coefficient exactly as was 
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the induced drag due to basic loading. Values of loading 
at the four spanwise points are found for the particular 
value of Cl and substituted in the appropriate form of equa- 
tion (9). If the more usual expression for induced drag as 

G 

a function of Cl is desired, then values of — 5 are substituted 


in the following expression which differs only algebraically 
from equation (9) : 


Cl 2 ~2c l<x 2 L\ « / 

— ^ (Y.056 — + 0.789 — ) 
a \ a a / 


fo.733 — + 0.845 — ) 

a \ a a /_ 


(13) 


Spanwise center of pressure. — Substitution of the values 
for individual loads in equation (Al7) of the appendix gives 
an expression for the continuous distribution of additional 
load. Integration of the increments of bending moment 
about an axis lying in the plane of symmetry will give total 
root bending moment. Then, with knowledge of the total 
load, an expression giving the spanwise location of the .center 
of load (on the quarter-chord line) can be found. In terms 
of the values of the dimensionless circulation, this expression 
is (from reference 5) 


Vc. 


G, Go G 3 (? 4 

0.352 -2+0.503 — 5 +0.344 — 2+0.041 — ? 

a a a a 

v ' « G, GZ Gs ~G~< 

0.383 -2+0.707 —2+0.924 —2+0.500 —2 


(14) 


Either of the two loading coefficients can be used directly 

G 

in place of the values of —• 


Pitching moment. — The value of the pitching moment due 
to additional-type load can be found exactly as was that due 
to basic load. However, account must be taken of the fact 
that the value of the moment is directly dependent on both 
the position of the reference axis and the amount of addi- 
tional load. This can be accomplished directly by using 
the results of the previous section which gave the spanwise 
position of the center of pressure. Since the center of pres- 
sure is placed by the method on the wing quarter-chord line, 
it is simply a matter of geometry to locate the longitudinal 
position of the center of pressure with respect to any reference 


axis. 


The value of the moment-curve slope 



is then 


simply the difference between the longitudinal position of this 
reference axis and the longitudinal position of the center of 
pressure expressed as a fraction of the mean aerodynamic 
chord. 


PROCEDURE FOR DETERMINATION OF AERODYNAMIC CHARACTERISTICS 
DUE TO GROSS LOADING 

Gross loading must be considered when a twisted wing is 
carrying other than zero net lift. Since the load distribu- 
tion is the sum of a constant and a variable, its shape varies 


with lift coefficient and thus the resulting characteristics 
must generally be determined at each lift coefficient. For 
the most part, excepting possibly the effect on the character 
of wing stall, the effects of twist on aerodynamic character- 
istics are of minor importance compared to the effects of 
sweep, aspect ratio, or taper ratio. However, where the 
twist varies with load due to the elastic properties of the 
wing, careful consideration must be given the gross loading 
if a true picture of the characteristics is to be obtained. 

Span load distribution. — Since basic and additional load 
are additive, the gross load distribution at a given C L is 
simply the sum of the two. The magnitude of the total 
load is equal to the value of the additional load. Thus, the 
gross load coefficient is given by G„ = G no -\-G na . 

It is sometimes of importance to determine the twist 
required in a wing to have a given gross load distribution 
at a particular lift coefficient. If this distribution is chosen 
and the additional load distribution for the wing correspond- 
ing to the particular lift coefficient is subtracted from it, 
then the necessary basic load distribution is defined. As 
noted previously, the twist required is readily determined. 
For wings having pronounced sweep or very low aspect 
ratio, the twist required to give load distribution approaching 
the desired elliptical may become great even where the 
distribution is desired at a lift coefficient of only 0.2 or 0.3. 

Lift-curve slope. — If the twist is constant, the lift-curve 
slope for the wing is that found for the untwisted wing. 
However, if the twist varies with load, then the lift coefficient 
must be determined for each angle of attack in order to fi n d 
the lift-curve slope. The required values of the lift coeffi- 
cient can be determined from equation (3) once the values 
of G n are known for each angle of attack. 

Induced drag. — Total induced drag can be found by 
substituting the values of total loading coefficient in equa- 
tion (9) . It is, of course, impossible to make C Di a function 
of C L when twist is present since C D is a function of load 
distribution and this varies with each lift coefficient. (Note: 
Cd^Cd -\-C Di .) 

7 *0 l a 

Pitching moment. — The gross pitching moment is found by 
adding directly the pitching moment due to basic lift to the 
pitching moment due to additional lift. This gross char- 
acteristic must be carefully examined if twist is a function of 
total load or dynamic pressure. In these cases a new twist, 
and the pitching moment due to basic lift must be found for 
each flight condition. It is possible the resultant changing 
value of basic load moment can be of sufficient magnitude 
to seriously affect the stability of an airplane. 

CORRECTION FOR SECTION LIFT-CURVE SLOPE 

The development of the basic theory involves the assump- 
tion that each section on the wing maintains a lift-curve 
slope of 2 ir. This is implicit in the choice of the three- 
quarter-chord line as the location for the control points. It 
may be, desirable at times, however, to consider the effects 
of deviation in lift-curve slope from the theoretical value of 
2x. How this can be done by moving the control point from 
the three-quarter-chord line is shown in the following 
discussion. 
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The logic underlying the choice of the three-quarter-chord 
line as a location for the control points can be shown as 
follows. The velocity w induced at any distance h from an 
infinite vortex of strength I’ can be expressed as 

W= irh 


—The circulation is in turn related to the section lift coefficient 
c t of a section through the expression 

l= P TV=Ci\pV 2 c 


From these two expressions it can be found that 

-"(fX?) 

It can be seen that the term y indicates a change in the 

direction of flow with respect to the free-stream direction, 
which change decreases with distance from the lifting line. 
If the lifting line is assumed to be replacing a plate insofar 
as lift is concerned, then simplified lifting-surface theory 
requires that, at some distance from the lifting line, the 
direction of flow must be parallel to' the plate. In. effect, 
then, the induced downwash angle becomes equal to the 
angle of attack of the plate. It remains to determine how 
far aft of the lifting line the downwash angle must be meas- 
ured in order to properly relate the increase in circulation 
to the change in angle of attack of the plate. From the 
foregoing expressions and assuming small angles, the follow- 
ing relation can be written: 



Now, if a section lift-curve slope of 2ir is assumed, it is 
evident that must equal ^ or that h is equal to ^ c, thus 
fixing the point where no flow occurs through the plate at 
7; c aft of the lifting line or on the three-quarter-chord line. 


It follows directly that, if a lift-curve slope is less than 
2 x, then the control point should move forward of the three- 
quarter-chord line, and if greater, aft. Thus, adjustment 
can be made for changes in section lift-curve slope at each 
of the four control stations. The procedure for doing this 
will be shown in the application of the method. It is not 
clear just how well this procedure will account for changes 
in section lift-curve slope resulting from separation of flow. 
Therefore, some caution should be used in interpreting 
results obtained from this method where large changes in 
lift-curve slope from the theoretical are involved. 

In using this lift-curve-slope correction, it must be real- 
ized that an attempt is being-made to impose considerations 
of section characteristics on a theory which cannot rigor- 
ously allow such considerations. The question of the ap- 
plicability of the correction becomes of particular import- 
ance where the wing has large values of sweep of the quarter- 
chord line. As the subject theory is developed, the sec- 


tions containing the control point and the point at which 
load is determined lie in a line parallel to the free-stream 
direction. Thus any considerations of section character- 
istics can be applied only to sections lying along these lines. 
However, simple sweep theory indicates that it is more 
correct to make the controlling airfoil section that one lying 
along a line normal to the wing quarter-chord line. These 
two sections are quite different, of course, if sweep is present. 
If no taper is present, their thickness distributions are the 
same, but their thickness ratios vary by the cosine of the 
angle of sweep. If taper is present, both thickness distribu- 
tion and thickness ratio vary. 

Study of the problem has so far indicated that neither the 
approach of simple sweep theory nor of the subject theory 
is correct insofar as choosing the controlling airfoil section, 
but rather that the controlling section lies somewhere be- 
tween these two limits and varies with span position. In 
view of this, therefore, it should not be expected that true 
account can be taken of changes in airfoil-section character- 
istics. However, it is believed that use of the procedure 
proposed herein will indicate the trend of the changes in 
over-all wing characteristics to be expected from changes in 
section characteristics. Until such time as more detailed 
analytic and experimental studies of the problem are com- 
pleted, it is recommended that the results of applying this 
correction be considered largely qualitative. 


EFFECTS OF COMPRESSIBILITY 

The Prandtl-Glauert rule, which accounts for the effects 
of compressibility, is directly applicable to the subject 
simplified lifting-surface theory. The approximations and 
limitations of the Prandtl-Glauert rule are well known and, 
hence, no discussion of them is given herein. However, for 
convenient reference, the basis for correcting the predicted 
span loadings and the theoretical relation of these correc- 
tions to the simplified lifting-surface theory is presented. 6 

The Prandtl-Glauert rule simply states that, as the Mach 
number is increased, the span load distribution of a wing 
distorts as though the x dimensions of the wing were in- 
creasing as the ratio of Thus, the effects of Mach number 

on a given wing can readily be considered by finding the 
span loading at zero Mach number of a properly distorted 
wing. It can be seen that increasing the x dimensions of the 
wing results in increasing the angle of sweep and increasing 
each local chord (or, in effect, a decrease in aspect ratio) 
while leaving the span and taper unaffected. Thus 




A^=tan _1 


tan A 
6 


The foregoing is applicable to any type of span-loading 
theory. The distorted wing is considered for the Mach 
number in question and the span load distribution found. 
This load distribution is then considered as being carried 
by the undistorted wing in order to find those characteristics 
dependent on span loading. 

* For the limiting case of Mach number equal to 1, the reader is referred to NAC A TN 1824* 
“Linearized Compressible-Flow Theory for Sonic Flight Speeds,” by Max. A. Heaslet, 
Harvard Lomax, and John R. Spreiter. 



8 


REPORT NO. 921 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


In the case of the subject span-loading theory, the compressi- 
bility correction can be considered in another way. As before, 
the ang le of sweep is effectively increased by increasing Mach 
number. However, the increase in local chord in the 

ratio i can be regarded as an increase in the distance between 

the lifting vortex and the control point on the undistorted 
wing. From the previous section it will be recalled this 
corresponds to an increase in section lift-curve slope — in this 

case exactly in the ratio of Thus, the theoretical section 

lift-curve slope, where compressible effects are included and 
the control point is held at the three-quarter-chord point, 

implicity becomes ^ Any corrections for section lift- 

curve-slope change thus become based on the deviations of 
the experimental section lift-curve slope at a given Mach 

2jt 

number from the value of -q at the same Mach number. 

It can be seen that this increase in c !a is exactly equivalent 
in effect to the decrease in aspect ratio. It is obvious from 
this that two techniques can be used in applying the com- 
pressibility correction to the simplified lifting-surface theory. 
Either the wing geometry can be appropriately altered or 
the disposition of the original layout of the lifting vortex 
and control points can be altered. The choice of procedure 
is governed entirely by the relative simplicity in handling 
the computations. As will be shown under the following 
sections on the application of the method, both are used, 
the choice depending on the type of loading and wing being 
considered. 

APPLICATION OF METHOD 


LOADING CHARACTERISTICS FOR ARBITRARY WINGS 


also possible. It can be shown (see appendix) that, if the 
angle of sweep is fixed, all the influence coefficients become a 

function of — alone. Thus, if a complete set of coefficients 

Cy 

for one taper ratio is determined throughout the aspect-ratio 
range, the coefficients for all other plan forms having the 
same sweep can be related to these coefficients by relating 

the values of — of the wing in question to those for which the 

Cy 

coefficients were determined. Further consideration shows 
that the effect of section lift-curve-slope deviation from the 
theoretical values can also be accounted for by a change in 

the value of — in the ratio of experimental to theoretical 

Cy 


section lift-curve slopes. Finally, if the compressible-sweep 
parameter ^A^=tan -1 is used rather than the geo- 

metric sweep angle and the value of — is adjusted by the 


factor p, then the effective plan-form change due to com- 
pressibility is accounted for. 

To simplify the use of the method, therefore, the influence 
coefficients for symmetric loading a, n for the seven-point 
method have been computed and are presented in table I 
and figure 1 as functions of the compressible sweep parameter 
and the parameter H v . For most wings the values of the 
coefficient a„„ can be obtained directly from figure 1. How- 
ever, for certain wings of extreme plan forms the values of 
H v will be such that the values of a„„ will lie off the charts. 
Where linear extrapolation of the curves is not acceptable, 
table I provides sufficient values of a,* to enable extension 
of the curves. Adequate accuracy of the final results will 
be obtained if the value of a vn is read to two decimal places. 

The parameter H v is defined as follows: 


In the foregoing sections, a method for predicting the span 
loading on wings has been outlined and it has been shown 
how other characteristics can be found from the span loading. 
Further, it has been shown how corrections can be applied 
to approximately account for deviations of the section lift- 
curve slopes from the theoretical value of 2 it and for the effects 
of compressibility. 

With the information thus far given, it is possible to predict 
the span loading and resulting characteristics of an arbitrary 
wing from a knowledge of the wing geometry only. Appli- 
cation of the procedure, however, shows that the most time- 
consuming portion is that of computing the influence coeffi- 
cients a, n , sixteen of which are required for the seven-point 
method. Examination of the theory shows that, if the num- 
ber of control stations is fixed, the influence coefficients 
become a function of the wing geometry only, that is, sweep 
and chord distribution. It becomes immediately apparent 
that, if the number of control stations is chosen, then the 
corresponding influence coefficients can be presented in 
graphical or tabular form as a function of wing geometry. 
Thus, the greater part of the computing work associated 
with the method can be eliminated, since the same coeffi- 
cients are used to find any form of symmetric loading for a 
given plan form. 

Further simplifications, not so immediately evident, are 


Hy = dy~ 

Ky 


_b 

Cy/P 


where 

d v scale factor which is given on each chart 


Cl 


K ' 2tt/p’ 


the ratio of the two-dimentional experimental lift- 


curve slope for the airfoil section at station v to the theo- 
retical value, both for the Mach number under con- 
sideration 
P VI— M 2 

— ratio of wing span to the chord of the wing at the span- 

Cy 

wise station corresponding to the control point v 
The application of the method to the case of the arbitrary 
wing plan forms can be outlined as follows: 

1. Compute the value of H r at each spanwise station v 
and the value of A3, both for the Mach number in question. 
(Note that the effects of section lift-curve-slope change and 
of Mach number are completely accounted for within the 
limits of the method by these parameters, and the predicted 
results will include these effects.) 

2 . From figure 1 or table I and with the value of H v , find 
the values of each of the 16 <z„„ coefficients. 

3 . Insert the values of a, n in equation (6) for basic load- 
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ing or equation (11) for additional loading and solve the 

simultaneous equations 8 for the values of 0 n or ^ n 5. 

a 

G 

4. With the values of G„ or known, other wing char- 
es 

acteristics can be found by substituting these values in the 
appropriate equation as indicated by the following table : 


Basic loading Additional loading 

C D{ Equation (9) Equation (9). 

C D JC L 2 Equation (13). 

Vc.p. Equation (14). 

C m Equation( 10) 


ADDITIONAL LOADING CHARACTERISTICS OF STRAIGHT-TAPERED 
WINGS 

The previous section has shown how to apply the method 
to determine aerodynamic characteristics of an arbitrary 
wing. For the more common case of straight-tapered wings, 
it is possible to prepare charts giving wing aerodynamic 
characteristics due to additional loading directly as a func- 
tion of wing sweep, aspect ratio, and taper ratio. 

The method has, therefore, been used to find the additional 
span loading and some of the wing characteristics for a series 
of wings of varying plan forms and having constant section 
lift-curve slope across the span. Some 200 wings (the range 
of plan forms, but not the total number of wings is shown in 
fig. 2) were included in the study, the results of which are 
presented in figures 3 to 6. Each of these figures is a chart 
giving the variation of a particular wing characteristic with 
wing sweep and aspect ratio for five values of taper ratio. 
Figure 3 gives the value of the local loading coefficient at 
each of the four spanwise stations. Figure 4 gives the 
value of wing lift-curve slope. Figure 5 gives the spanwise 
location of the center of pressure. Figure 6 gives the chord- 
wise location of the center of pressure or the aerodynamic 
center 7 measured from the leading edge of the mean aero- 
dynamic chord. 

To find the desired characteristic for a given wing, the 
chart for the proper taper ratio is entered with the com- 
pressible sweep parameter (^Ap= tan -1 and the de- 

sired value read from the curve- — -or interpolated — for the 

a A 

proper value of the aspect-ratio parameter Thus, it is 

possible to find, with no computation, many of the charac- 
teristics for untwisted wings. It should be remembered 
that, since basic and additional loading are considered 
independent, these characteristics not only represent the 
gross characteristics of untwisted wings, but also .the addi- 
tional loading characteristics of twisted wings. They are 
thus applicable, to some extent, to all wings having straight- 
taper and constant section lift-curve slope across the span. 


• Experience has shown that the method given in reference 6 for the solution of a set of simul- 
taneous equations is most satisfactory. 

7 For the case of the straight-tapered wing, this chordwise location of the aerodynamic 
renter can be simply expressed in terms of the geometry of the wing and the spanwise location 
of the center of pressure. Thus 


905385—51 39 


“• '-i+sf™ I>-’-.-3TO] A tan A 


TABLE I.— SYMMETRIC INFLUENCE COPIFFICIENTS, a,„. 
WHICH LEAD OFF THE CHARTS OF FIGURE 1 
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DISCUSSION 

COMPARISON WITH OTHER THEORY AND WITH EXPERIMENT 

The extent to which the method can be evaluated varies 
widely with the particular characteristic in question. In 
general, the accuracy of predicting the characteristics asso- 
ciated with basic loading can be evaluated in only a limited 
manner, since few swept or low-aspect-ratio wings having 
appreciable twist or camber have been experimentally or 
theoretically studied. On the other hand, many wing plan 
forms have been studied and from those results it is possible 
to assess fairly well the accuracy of the method with regard 
to additional loading. 

Basic loading. — The prediction of the effect of twist and/or 
camber on load distribution can be evaluated to a degree 
by comparison with other theories. Such a comparison is 
given in figure 7 wherein the predicted loading given by the 
subject method is compared with those given by the method 
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developed by V. M. Falkner and that of reference 7. It may 
be assumed tbat the Falkner method is the most accurate of 
the three since it approaches a true lifting-surface theory. 
From the comparisons shown in figure 7, it would appear 
the subject theory is nearly as good as the more complex 
theory and considerably better than the lifting-line theory 
of reference 7. It must be recognized that twist introduces 
lar ge induction effects even on high-aspect-ratio wings. 
Therefore, it follows that the deviations of the predictions 
of the subject theory and those of the theory of reference 7 
from those of Falkner’s method are indicative of the ability 
of the two methods to account for the induction effects. 

A check on the ability to predict Cm 0 and a Tf) can be had 
from comparison with the experimental results shown in 
reference 7. Here a wing having 30° of sweepback was 
tested without twist (wing designated 24-30-0) and with 
8.5° of twist (wing designated 24-30-8.5). The force tests 
show that the twist shifts the angle of zero lift 2.6° while 
the present theory gives a value of 3.3°. The force tests 
indicate that twist produces a value of C m(i of 0.05 which 

is the same as that predicted by the theory. 

A further comparison of experiment and the present 
theory can be made using results of pressure-distribution 
tests of a wing having 60.8° of sweepback. The geometrical 
characteristics of the wing together with the experimental 
results and theoretical predictions are shown in figure 8. 
The comparisons are made for zero lift where only basic 
loading exists. It can be seen that the agreement is good 
even for this unconventional plan form. 

On the basis of the foregoing comparisons, it seems justifi- 
able to conclude that the subject method can adequately 
predict the effects of twist and/or camber on the character- 
istics of wings of arbitary plan form. 

Additional loading.— Much experimental data and several 
theories are available for comparison with the subject theory 
in regard to the prediction of additional-type lift. Figure 
9 shows a comparison of the predicted loadings and those 
obtained experimentally for six wings varying in sweep, 
aspect ratio, and taper ratio. Figure 10 compares the 
variation of lift-curve slope with aspect ratio as predicted 
by the subject theory, by two more rigorous theories (refer- 
ences 8 and 9) , and by a theory directed at the limiting case 
of zero aspect ratio (reference 10). Figure 11 compares 
experimental and theoretically predicted lift-curve slopes 
(assuming a section lift-curve slope of 2jr per radian) for 
two families of plan forms covering a wide range of aspect 
ratios (references 11 and 12). The comparison between 
experimental (references 12 and 13) and predicted results 
for a random group of wing plan forms is shown in figure 
12 for lift-curve slope and figure 13 for aerodynamic-center 
location. All experimental values of C La and aerodynamic- 
center location were measured at zero lift. It is evident 
that in almost every case the method gives an excellent 
prediction. 


From all of this comparative material it seems possible to 
safely draw the conclusion that the subject method can 
satisfactorily predict the additional-type span loading on 
wings of arbitrary plan form. Further, it seems possible 
that equally good predictions can be made of those charac- 
teristics primarily dependent upon such span loading. 


EFFECT OF PLAN-FORM VARIATION ON THE BASIC-LOADING 
CHARACTERISTICS 


To study the effects of plan-form variation on the basic- 
loading characteristics of linearly twisted, straight- 
tapered wings, the characteristics of a representative group 
of wings (see shaded wings, fig. 2) having unit washout 8 
have been computed and are presented in figures 14 to 19. 
The Mach number was taken as zero and the section lift- 
curve slope as 2ir. The basic-loading characteristics con- 


C i-C 

sidered are the loading coefficient -2- (figs. 14 to 17), the 

€C a v 

c m 

pitching-moment coefficient due to twist (fig. 18), and 


the angle of attack of the root section for zero net lift — 
(fig. 19). 

Magnitude and spanwise distribution of load.' — Examina- 
tion of figures 14 to 17 reveals that the aspect ratio influences 
only the magnitude of the load and is in fact the predominate 
influence on the magnitude of the load. Reductions in as- 
pect ratio from 6.0 to 3.5 and 1.5 result in approximately 
35-percent and 70-percent reductions, respectively, in load 
due to twist for either the unswept or 45° swept-back wings 
(fig. 15). 

Sweep, either forward or back, tends to reduce the magni- 
tude of loading, although appreciable reductions are produced 
only by sweep angles greater than 45° (fig. 14). Sweep also 
affects the load distribution such that the load on the outer 
section of the wing is shifted inboard by sweepforward and 
toward the tip by sweepback; as will be seen, this is similar 
to the effect of sweep on the additional-type load'ng. Since 
an increase in aspect ratio magnifies the loading, it also mag- 
nifies the effects of sweep on the loading as is shown in 
figure 14. 

As shown in figures 16 and 17, variation in taper ratio has 
little effect on the magnitude of basic loading; and variations 
in taper ratio, for taper ratios larger than 0.5, have little 
effect on the load distribution. However, for taper ratios 
less than 0.5, the loading on the outer section of the wing 
shifts inboard. These effects of taper ratio on loading are 
magnified by increases in aspect ratio. 

Pitching moment. — That the pitching moment due to 
twist is primarily a function of sweep and aspect ratio is 
shown in figure 18. The magnitude of the pitching moment 
increases as either aspect ratio or sweep is increased so that 
pitching-moment coefficients as large as 0.008, for 1° of twist, 


* In this case, 1° was chosen, and for any larger amount of twist the effects are proportional. 
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exist on wings having large aspect ratios and sweep angles. 
The effect of taper ratio is relatively small, the greatest 
being evidenced at the small values of taper ratio. For 
example, reducing the taper ratio from 0.5 to 0 reduces the 
pitching moment due to twist about 30 percent. 

Angle of zero lift. — Although the effects of plan form on 
the angle of zero lift a ro may not be very important, some of 
the trends indicated in figure 19 are of interest. For the 
range of plan forms' represented in figure 19, the angle of 
zero lift varied only about 20 percent from the mean value. 
This is small compared to the effects of plan form on the 
magnitude and distribution of loading and on the pitching 
moment. In contrast to the small effect of taper noted 
previously, taper ratio appears to be the predominant influ- 
ence on a r0 particularly at large aspect ratios and large 
sweepback, The effect of aspect ratio and sweep are sec- 
ondary but not negligible. 


EFFECT OF PLAN-FORM VARIATION ON THE ADDITIONAL LOADING 
CHARACTERISTICS 

Span load distribution. — To show directly the effects of 
wing plan form on the additional span loading, figure 20 
has been prepared by cross-plotting the data of figure 3. 
It shows that increasing the angle of sweepback or the taper 

ratio serves to move the loading ^as defined by ou ^' 

board. 

It can be noted in figure 5 that the spanwise center of 
pressure is independent of aspect ratio for certain combina- 
tions of taper ratio and sweep angle. These values of taper 
ratio are plotted against sweep angle in figure 21. Further, 
for the wing geometry represented by the curve of figure 21, 
the loading is approximately independent of aspect ratio 
(fig. 3), and also is approximately elliptical. For elliptical 
loading, the loading coefficient at the four spanwise stations 
have the values 1.273, 1.176, 0.900, and 0.487. It can be 
seen that these values compare closely with those given in 
figure 3 for the plan forms specified in figure 21. The farther 
the wing geometry departs from the configuration represented 
by the curve in figure 21, the greater the change of loading 
with aspect ratio and the rate of distortion from an elliptical 
load distribution. Also, examination of figure 3 will show 
that elliptical loading cannot be maintained, when a wing is 
swept, by altering aspect ratio alone; however, all wings 
approach an elliptical loading as aspect ratio approaches zero. 

The wings specified by the curve of figure 21 have the 
property that their aerodynamic characteristics can be 
expressed in a simple manner similar to the case of unswept 
wings with elliptic plan forms. Namely, the induced drag 

q d i 

is given approximately by spanwise center of 

4 

pressure by and aerodynamic-center location with 


respect to the mean aerodynamic chord by 


1 i 0.342 — 0.567X — 0.908X 2 
a ' C ' 4 -1 " 10(1 + X+X 2 ) 


A tan A 


From figure 4, it can be shown that the wings specified by 
the curve of figure 21 give the maximum lift-curve slope for 
a given aspect ratio and sweep angle. For example, see 
figure 22 for the case of A=3.0. 

In summary, the curve of figure 21 defines wings having 
the following approximate characteristics-: span loading 
distribution independent of aspect ratio, minimum induced 
drag, and maximum lift-curve slope for a given sweep and 
aspect ratio, a constant spanwise center of pressure, and an 
aerodynamic-center location that is a simple function of 
wing geometry only. 

Lift-curve slope. — Examination of figure 4 reveals certain 
general trends with respect to the effect of wing plan form 
on wing-lift-curve slope. For wings of high aspect ratio, 
the angle of sweep has a marked effect on lift-curve slope, 
with the maximum effect occurring for wings of infinite 
aspect ratio when the lift-curve slope is directly a function 
of the angle of sweep. As the aspect ratio approaches very 
low values, the lift-curve slope for the unswept wing is 
greatly reduced and the effects of sweep become small except 
for very large angles of sweep. Also, it can be seen that at 
very large angles of sweep the effects of aspect-ratio varia- 
tion on lift-curve slope become small. 

To better illustrate the separate effects of aspect ratio 
and taper ratio, the data from figure 4 have been cross- 
plotted to show the variation with aspect ratio of the lift- 
curve slope for various values of taper ratio and sweep angle. 
These results are shown in figure 23. This figure shows 
clearly how increasing the angle of sweep decreases the varia- 
tion of lift-curve slope with aspect ratio. It shows further 
that while taper ratio as compared to aspect ratio has only 
a small effect on the lift-curve slope of an unswept wing, 
taper ratio has a predominant effect on the lift-curve slopes 
of highly swept wings of moderate to high aspect ratios. 
For very small aspect ratios (i. e., 1.5), however, the lift- 

curve slopes of all the wings converge and become almost a 
linear function of aspect ratio, being essentially independent 
of the effects of sweep and taper (reference 10). 

Aerodynamic center. — In figure 6, variations of aerodyna- 
mic-center location from 15-percent mean aerodynamic 
chord to 45-percent mean aerodynamic chord are indicated 
for the range of plan forms studied. It will be noted that, 
for taper ratio X=0, the aerodynamic center moves aft for 
sweepback and forward for sweepforward. At taper ratios 
of 1.0 and 1.5 the aerodynamic center moves forward for 
sweepback and aft for sweepforward. For X=0 the effects 
of aspect ratio are largely confined to the swept-back wings 
and as taper ratio is increased the effect of aspect ratio 
decreases for swept-back wings and increases for swept- 
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forward wings. Where the effects of aspect ratio are signifi- 
cant, an increase in aspect ratio generally moves the aero- 
dynamic center aft. 

CONCLUDING REMARKS 

The simplified lifting-surface theory presented herein 
enables the rapid, accurate prediction of symmetric span 
load distribution for wings which have symmetry about the 
root chord and a straight quarter-chord line over the semi- 
span tlnd which can have arbitrary chord distribution, 
sweep, aspect ratio, and continuous twist. Modifications 
to the method are shown by means of which approximate 
account can be taken of changes in section lift-curve slope 
and by means of which the effects of compressibility as 
predicted by the Prandtl-Glauert rule can be included. 


With the charts presented in the report, the load distribu- 
tion can be obtained directly for many wings and can be 
obtained with slight additional computation for all wings 
falling within the limitations prescribed. 

From the span loadings determined by the method, it is 
shown how several important wing characteristics can be 
determined. Comparison With experimental results in- 
dicates that the wing characteristics so obtained are reliable 
for a wide range of plan forms. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., 194 - 7 - 1948 . 


APPENDIX 


The mathematical process of finding the loading distribu- 
tion for wings evolves itself into two parts: first, the deter- 
mination of the integral equation which relates the downwash 
at a given point to the integrated effects of the bound vortex 
and trailing vortex sheet; and, second, the solution of the 
integral equation to determine the unknown loading distri- 
bution factor. 


ADDITIONAL SYMBOLS USED IN THE APPENDIX 




b v y 

byn 

By 


! 


Fourier coefficients of the loading distribution series 


mathematical series coefficients 


&nk 

fnii \ 

In „ ) 

Qvn ^ 
Qvv 
Qvn ^ 

k 


L(v, v) 



M 


r 

s 

x 

X 


y 

V 

«i 

V 


Bi,8 a 




loading interpolation-factors 
mathematical series coefficients 


downwash influence coefficients 


value pertaining to interpolated span station 

G=cosy) 


| downwash integrand functions 


number of span stations taken to numerically 
integrate the downwash integrand function 
absolute distance from the downwash point to an 
elemental vortex, feet 
distance along the load vortex, feet 
longitudinal coordinate, positive downstream, feet 
longitudinal coordinate pertaining to the load 
vortex, feet 

lateral coordinate, positive to the right, feet 
lateral coordinate pertaining to the load vortex, 
feet 

induced angle in the wake of the wing, radians 
dimensionless lateral coordinate pertaining to the 


load vortex 



positive angles between a vortex line and the lines 
joining the ends of the vortex line and a down- 
wash point, radians 


integer pertaining to 


span station y ^17= cos 


Hi integer sequence of the Fourier series for loading 

distribution, also pertaining to span station n\ 

G= cos ^¥i) 

DETERMINATION OF THE DOWNWASH INTEGRAL EQUATION FOR SWEPT 

WINGS 

Downwash induced by the trailing vortex sheet. — 



The downwash due to the trailing sheet is, for an arbitrary 
elemental trailing vortex, given in several references (e. g., 
see Glauert, reference 14). The downwash at a point in the 
xy plane is given by 

(cos 0i + cos 0 2 ) (Al) 


where the induced velocity is positive for downwash, 0 t and 
0 2 are the inside angles between the trailing vortex line and 
the lines from the ends of the trailing vortex line to the 
downwash point, h is the perpendicular distance from the 
downwash point to the trailing vortex line, and dr is the 
strength at a given span station of the trailing vortex. From 
the foregoing sketch it can be seen that 

h=y—y 

cos 01 — V(z— 5) 2 +(y— y) 2 


cos 0 2 =1 
and 


(since 0 2 =O with the trailing vortex extended to 
infinity) 

x= |y | tan A 


Substituting these values into equation (Al) and integrating 
will give w„ due to trailing sheet equal to 


J_ rw 1 I" x— jyltanA 

4 r J-6/2 y—y L 1 + V (*— \V\ tan A)*+ (y—y) 
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Downwash induced by the load vortex.- 


0 Vcos 2 7v +2 tan A ~ y ) T '^ dy 
- 6/2 (x+y tan A) V (x+y tan A) 2 + (y—y) 2 



ds r- -Bound or 
/ load vortex 


Trailing 

vortex 

sheet 


The downwash generated by the swept load vortex is the 
sum of the downwash generated from each semispan. Along 
the load vortex the circulation is continuously varying. The 
downwash at (x, y ) due to the small element ds (see accom- 
panying sketch) is given by (e. g., see Glauert, reference 14) 

(A3) 

For a swept load vortex, for the right semispan 


h=x cos A — y sin A 
r= \ (x—x) 2 +(y—y) 2 
x—'\y\ tan A 

while for the left semispan 

h= x cos A 4 -y sin A 

and s, r, and x are the same as for the right semispan. 
Substituting these values into equation (A3) and integrating 
will give w xy due to load vortex equal to 

J_ fo (x cos A + 7 / sin A) I' (y)dy , 

44 J -t /2 cos A [(r— T/[ tan A)-+ {y— y) 2 Y 12 

1 Cm (; x co s A — ;/_s i n A) F (y)dy .... 

4ir Jo cos A[(x— jyj tan A) 2 +(j/— y) 2 ] :i/2 

Total downwash.-— The total downwash due to the trail- 
ing sheet and the load vortex is equal to the sum of equa- 
tions (A2) and (A4), or 

1 C h '' 2 1 r . r— lyj tan A 1 , 

w„= r — = 1 r 0/) f ty+ 

4ttJ - i ,/2 y— y L \{x— |y| tan A)-+(j/— y) J 

l_ f° (x+y tai l A) V ( y)d y 

4Tj-m [(x—'lyi tan A) 2 + (y— f)' 2 f /2 

1 P /2 (r— y ta n \)T(y )dy . 

4ttJo {(x-\y\ tan A++(y— y) 2 ] m 

The last two integrals of equation (Ao) can be integrated 
by parts and put in the form of the first integral. With 
/ 

± 2^=0, equation (A5) becomes 
■I- ( C bl2 1 fi i x—\y\ tan A "1 n r\ j- 

Wxy 4tt| J-6./2 ?/— yL |y| tail A)-+ (?/ — y) 2 J V V 


.cos 2 A 


-x tan A 


-y) C'Cy)dy 


Jo {x—y tan A) V (z— V tan A) 2 + (y~~y) 2 

rn 

Equation (A6) can be simplified algebraically and rewritten 
into the form 


1 r* 

4xJ-6/2 y—y 

1 C° V(^+y tan A y +(y-w r ^ d - . 

4x J —bj 2 (x+y tan A) (y—y) V t 

1 rw Vfr-y tan S+ 5EE r@#+ 

4tt Jo (x—y tan A) (y—y) a 

2 tan A V^+y 2 .. . 

4x(x 2 -y 2 tan 2 A) 1 [ ) { ) 

The three integrands of equation (A7) have an infinite 
point at y—y. The integrand in the second and third inte- 
grals can be made continuous by subtracting the function 

r'(v) " r'(w) 

-■ Equation (A7) remains balanced if =■ is also added 

y—y v-y 

to the first integrand. As will be seen presently, the first 
integral, the integrand of which retains the infinite point, will 
give a finite value after integration. Equation (A7) becomes 
(adding the fourth term to the second integral) 

2_ F 6/2 C'(y)dy _1_ j”° [~ V(s+y tan A) 2 +(y— y) 2 
Wxl 4x J - 6/2 y—y ' 47r J — 6/2 L (x + y tan A) (y y) 

_ 1 j. - tan A^ 2 +J/ 2_ j r , (n\ ( ^j+ 

y-T+ x’-y 2 tan 2 A J 1 

-}- r I (AS) 

4xJ« L (x—y tan A) (y—y) y—yj y 

Equation (AS) can give the downwash at any point in the 
xy plane. From the discussion of the theory in the text, the 
downwash must be found along the three-quarter-chord line 
of the wing. For a given span station the x coordinate of a 
point on the three-quarter-chord line is given by 

*=2 + |y| fan A 

The nondimensional equation for downwash angle at any 
span station on the three-quarter-chord line can be obtained 
by substituting this value of x into equation (A8) and by 
using the following nondimensional relations: 


__JL 

v ~b/2 


y 

"“6/2 
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for !J<0, 


& 72 =f+M tanA 



The resulting equation is 

(vH (A») 

where the function -£(17,17) is given by, 


T (n 17 ! — 1 ( V[l + (t/c) (\y\+v) tan A] 2 + (b/c) 2 (y—y ) 2 K 2 tan A y'[l (6/c) ln| tail A] 2 +(6/c)V 

V} ( b/c- ) {n-v) \ 1 + (b/c) (h| +17) tan A ) Ml + (b/c) (|,|-,) tan A][l + (b/c) (|,| +,) tan A] 

and, for rj > 0, 


L(v,v) 


1 ( V[l + (6/c)(M-i?) tan Aj 2 +(6/ C ) 2 h-h)* 
(b/c) ( v -v) \ 1 + (b/c) (h| -17) tan A 


j 


(A10) 


Equation (AlO) is complete for positive and negative values 
of the coordinate 17. In the case of the unswept wing where 
A equals zero, equation (AlO) simplifies considerably and 

the function £(17,17) is given by, 

for — l<j?<l. 


L(v,v) = (b/c)l n —q) [Vl + (^/c) 2 h— v ) 2 — lj (All) 

MATHEMATICAL SOLUTION OF EQUATION (A9) 

Arbitrary load distribution— The solution of the integral 
equation (A 9 ) to determine the unknown loading-distribu- 
tion function G(ij), was obtained by Weissinger (reference 1) 
using a method introduced and applied by Multhopp 
(reference 15 ) to the problem ’ of determining downwash 
at a straight bound vortex. Weissinger’s method consists 
of applying the boundary conditions that the flow shall be 
tangent to the plate at the three-quarter-chord line at a 
number of span stations, then performing a numerical 
integration of equation (A 9 ) at each station. The result 
is a set of simultaneous linear equations in which the 
unknowns are the values of G( r?) at the span stations 
chosen. 

Introducing in equation (A 9 ) the spanwise trigonometric 
variables </>, and <£ defined by 


and 


ij= cos </>, 
17= cos <j) 


gives 

/w\ _1 f * 


G'( 4 >)d 4 > b/c , i ' 17 

cos </>— cos <j>, 2ir Jo 


L(<f>v,<t>)G' (<t>)d<f> (A12) 


where b/c, is the value of b/c at span station v. 

In equation (A12) let G(</>) be given by the Fourier series 


as 


where 


G(<f>) = X) a. sin mi<£ 

e.=i _ 

I fax 

I G(</>) sin yi<f>d<t> 

1 TTJ 0 


(A 13 ) 


Since G(</>) is of interest only in the interval 0 to x ( — 6/2 < 
y<b/ 2 ) the» 


( Al4 > 


Multhopp in reference 15 develops a quadrature formula 
which is a simple analogy to Gauss’ (e. g., see reference 6) 
mechanical quadrature. This integration formula is a good 
choice for the functions represented by the series as given in 
equation (A 13 ) since it integrates exactly functions repre- 
sented by the trigonometric series to the 2 m th harmonic. 
The quadrature formula is given by 

J l ir m 

/(J7)di7=— S f(rin) sin </> n (A 15 ) 

“I A 71 = 1 

where </>= and J (r) n ) is the value of /(? 7) at 17— cos 
Equation (A 14 ) can be integrated by equation (A 15 ) to give 
2 m 

X G(<t>n) sin Ur<t>n (A 16 ) 

Substitution of equation (Al6) in (A 13 ) gives the new form 
of the loading-distribution function as 

2 771 

G(<f>) ~ ^rn X Gn S sin mi <t>n sin mi<£ (A 17 ) 

1 Mi=i 


where G n -=G(<f> n ). 

Then the derivative of equation (A 17 ) is given by 


dG(</>) 

d<t> 


G'( 4 >) 


2 ra 771 

— tt XI G n X! Mi sin /ii <p„ cos mi<£ . 

ni ~ x ~ 1 71=1 

(Ais) 


The first integral of equation (A12) containing the infinite 
point can be integrated directly with the series given by 
equation (A 18 ). Using the following integral derived in 
reference 14 

r cos n</> , _7r sin n<p 0 

Jo cos 4 >— cos </> 0 sin <j>o 

then 


1 C x G f (<f)d</> 

T Jo cos <f > — cos </>, 


2 

irfm+l) 


m m 

sin /u<£» 

n = l ju,=l 


r r cos mi 4 >d</> 

Jo cos </>-— cos </>, 


2 a mi sin m i<t>n sin mi<£, 

m-f 1 ^ sin <t> , 


m 

— 2 b„G , — X* 


n=l 


2 b, n Q n 


(A 19 ) 
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where the prime on the summation sign indicates that the 
value for n=v is not summed. Physically, equation (A19) 
gives the downwash angle at an infinite distance down- 
stream for any wing geometry. This downwash angle is 
equal to twice the downwash angle at the quarter-chord line 
for the unswept wing. 

The b„ coefficient of equation (A 19) are given by the 
following equation: 


For n=v 


' (m+ 1) sin 


^ ' /i i sm 


For n?±v 
b v 


Y (A 20 ) 


-1 


(m+1) sin^ftTi 


— Mi sin ni<f>n sin mi <+ 


Equation (A20) can be simplified by making the indicated 
summations. With the relation 

m / ™ \ 

X) Mi cos ^ 2 = real part of ( 2_, fitf’W ) 

jU,=l \Ml=l / 


1+m cos (m-j-l)a:— (m+1) cos mx 
2( — 1 + cos x) 


(A21) 


then for n=v 


and for n^v 

by n 


5 _ m+l 

v " 4 sin <j)y 

sin 4>n l-(-l)"-- 

(cos <j> n — cos <j> r y 2(m+l) 


Y (A22) 


• 

The integration of the second integral of equation (Al2) is 
done by a method similar to that used for the first integral. 
The integral of the product of the A (+,,<£) function and the 
terms of the trigonometric series of G' (</>) is too difficult for 
direct integration but can be done with the aid of the in- 
tegration formula. This integration can be made for an 
arbitrary number of points M which can differ from m. 

Using a trapezoidal integration formula for the trigono- 
metric variable gives 


/>***- ETI 0 


7 ( < t> o) 47 (4>m+ i) 


m n 

S/(+) 

/*= i J 


(A23) 


where <f>„— With the definition 
2 m 

S Ml Sin C ° S ^ ( A24 ) 

equation (A18) can be written in shortened form as 

m 

GJJy) (A25) 

71=1 

Applying the integration formulas of equations (A23) and 
(A25) to the second integral of equation (A12) gives 


G’ (4>)d<t>= £ QJn&)d* 

~Cyhi hn h{M+\)\_ 2 


M 

J A (< py } <p[i 
„=i 



(A26) 


Equation (A26) can be simplified as follows: Let/ n(1 =/„(<£,,) 

fXTT 

where /„(<#>„) is given by equation (A24) for <f> — 

and let 
Then define 


Summarizing the development for the computation of 
spanwise loading over the whole wing span, the loading at 
m span stations is found from the simultaneous solution of 
m equations, 


'2 (M+1) 


( 


2 


M \ 

Si vpfw ) 

e= l / 


(A27) 


A— j 

71 = 1 


Then equation (A26) can be written 

Jhi Lit,, 4>)G'(<t>)d4>=^ £ g vn G n (A28) 

Z7T JO to=1 

Equation (Al2) can now be put in a form which allows 
solution by simultaneous equations by substituting in it the 
equalities given by equations (A19) and (A28). Making this 
substitution equation (A12) becomes 

/nn\ m h m 

( ^) =2 byyGy-T,' 2 b, n G n ^ S 9*0* 

\ V Jr 71 = 1 71 = 1 

=^2 b y y-\-— g v ^ Gy — ^2* (^b yn ^ g 


where 


v=l, 2, 3, . . . m 


(A30) 



= angle of attack at span station v 


A yn =2byy+-~- g U y for n=v 
c v 


= — 2 b rn +~ g,n for n^v 

byv and b rn are from equation (A22) 
gw gvn for 71 v 


v=l,2,3, . . . m 


(A29) 


g vn is from equation (A27) 


where the prime on the summation sign indicates that the 
value for n—v should not be summed. 


b wing span 

c,~ chord at span station v 
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L,„ is obtained from equation (A 10 ) for the swept wing, and 
from equation (All) for the unswept wing, with 77 = cos 


771+1 


and 77 = cos and where /„„ is from equation (A24) with 


fiT 


M+ 1 

Symmetric load distribution. — The computations for the 
symmetrically loaded wing can be considerably reduced by 
altering the preceding equations and coefficients. For a 
symmetrically loaded wing the distribution of local angle of 
attack is symmetrical about the plane of symmetry, then 


OL v OL m J r \^. v \ 

G v z= v / 


(A31) 


With symmetric loading, only the odd values of mi in equa- 
tion (A13) contribute to the loading, then equation (A24) 
becomes 

2 A 


/, 


m + 1 Mi = l,3,5...odd 


Ml Sin fi+n cos 


(A32) 


Since the change of symmetric loading distribution with span 
is opposite in sign from one side of the span to the other, then 

f n,M+l — m and /* n, M+l 0 
2 

Furthermore, examination of equation (A32) indicates that 

J n(i == Jm + 1 — n,i i 

Then for symmetric loading equation (A27) becomes 


— 1 


M-l 
2 _ 


where 


9*n — 2 (M+ 1 ) fnii{L yil 
m+l 


(A 33) 


71=1 . . . 


or f .m+l 

Jw= Zfnn for — 2 ~ - 

—f 

— j n\i for 71 — 2 

—jfnn for and m=0 


for and g=0 


(A34) 


Equation (A32) gives the values of 

For symmetric loading, with the relations given by equa- 
tion (A31), equation (A29) becomes 


Wl+l 


(v) r{ a ~ + l *-) a -~ti *>■•) a • 

m+l 


*=1.2,3, . 

#05385—50 40 


(A3 5) 


(The prime on the summation sign indicates that the value 
for n— v is not summed.) 


where 


for 


=b. 


for n— 


m+l 


It should be noted that the equation for L, n simplifies some- 
what for the symmetrically loaded wing since, for symmetric 
loading, it is necessary to consider only positive values of V. 
Further, equation (A 10 ) can be written using only positive 
values of V through the following relation: 

L* L (ij J^j) A(t?, 7?) === I J V f l -Z+llf+1— M 

In summary, the foregoing analysis for the case of the 
symmetrically loaded wing gives 

Tn-f 1 
2 

OCy^^ ^_ 'i Q/vnGn 
71 = 1 


where 


r=l,2,3> • 


b _ 


m+l 


a yn — 2b yy +- g yy for n=v 

Cy 


= — 2 B yi >+^ g rn for n^v 


(A 3 6 ) 


(A37) 


g^=g yn for n=v 


M-l 
■ 1 2 


Syn 2 (M+l) 

/*m = 2 /„ m for 719^— 


=jnn for 71 = 


m+l 




=fn„ for m =0 

. Jnjx £ m+l n 

— 2 ^ — 2 ’ ^ ® 


™ I 1 +—J 

m-f- 1 Mi=I,3,5...odd 

m+l 
'4 sin <j>. 


X) Mi sin Hi<f>n cos Mitfv 


B yn =b yn +b,, m+l ~ n for 7i 


m+l 


=b yn for 7 i= 


m+l 


b 

L* 


(cos <t> 
1 


sin <f> H ri — ( — 1 )’ , ~ , ~] 

*— cos<£,) 3 |_ 2 (m+l) J 


"* b. 

-(m— m) 


I V^ 1 ^ tan A J* (m—v*)*— 1 1 
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— (Vr+Vn) 

V V 




1+2 — tan A ;. 

Cy 


2 tan A*y ( 1 +~ V? tan A 


1+2 — ,„ tan A 

Cy 




,, = COS ; 


, M = COS 


M+l 


~M + 1 


DEVELOPMENT OP THE PARAMETERS USED IN THE SIMULTANEOUS 
SOLUTION FOR SYMMETRIC LOADING 

The linear simultaneous equations of equation (A36) 

provide means for the solution for symmetric loading for any 
given symmetrical distribution of a y . The coefficients 

a yn , which are 2 ~ ^ in number, are functions of the wing 

geometric parameters of sweep angle A and the ratio — ■ The 

c 

principal work in a study of spanwise loading is to compute 
the a yn coefficients. For a study of a range of plan forms, 
these a rn coefficients can be plotted as functions of sweep angle 

and the ratio — ; however, — will vary spanwise due to the 

G C 

taper of a wing and plots of a vn as a function of — become 

unwieldy for a range of plan forms. A scale factor can be 

applied to — such that the spanwise variation will be effec- 
c 

tively nullified for a range of tapered wings. 

For a range of aspect ratio, the values of 7 for the out 

c 

board half of the wing semispan, ,>0.5, has maximum values 
for zero tapered wings (provided the plan-form edges are 
not concave) and for the inboard half of the wing semispan, 

77 + 0 . 5 , — has maximum values for the inverse tapered wings 
c • 

The ratio of — in the general case to these maximum values 
c 


of — provides a geometric parameter that varies approxi- 
mately as the aspect ratio. 

The inverse chord distribution for straight-tapered wings 
is given by 

k = yl ( 1 + X ) fAQON 

c 2[1 — h[(l— X)] 


For X=0 


For X=1.5 


b 1 

Ac 2 ( 1 -|,|) 

b 5 

Ac 2(2+|,|) 


The ratio of - to equations (A39) and (A40) gives, respec- 
tively, the new geometric parameter as 

Tyv — =2(1-,) ^ for 0.5 <tj<1 
\Ac />'=o 

(A41) 

•L_3S+iU ta0as , SM 

\dc/ l -1.5 


Let H be defined as 2/5 times the values of equation (A41). 
The factor 2/5 is introduced to give H the approximate 
values of a vn to simplify plotting procedures. Then equa- 
tion (A41) becomes, 

H= 4(i-77) i,o. 5£K r 

4(2+,) b [ (A42) 

~~25~ 

Introducing the effects of compressibility and section lift- 
curve slope as discussed in the text, equation (A42) becomes 
at the span station v, 


*-*©(&) 


where d, is a scale factor given by 

d for 0.5<,„<1 

= ~25*~ for 0<,,<0.5. 

For straight-tapered wings, equation (A43) simplifies to 
jj 2(1-,,) (1+X) /1\ . n 

H " 5 [ 1 -, A 1 - X )] U M,0 - 5 ^^< 1 

- (A44) 

Plots of a rn against H y in the range H„= 0 to 4, will give 
a yn coefficients for wings of any chord distribution with 
aspect ratios from 0 to 10 or 12 . 
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AERODYNAMIC CHARACTERISTICS FROM INTEGRATED SPANWISE 

LOADING 

lift coefficient. — The wing lift coefficient is given by 

Cl=A^ G(ri)dri=A-^ G(4>) sin 

or using the quadrature formula of equation (Al5), 

7 tA 


e L =- 


S Gn sin 4>n 


m+1 £=\ 

For symmetric loading equation (A45) simplifies to 


(A45) 


m— 1 


C ‘=STI (% + 2gA»m*.) (A46) 


the wing induced drag is given by considering the induced 
angle as one-half the value of the downwash at an infinite 
distance downstream. Then a t {<j > ) is given by one-half the 
value of equation (A 19). Equation (A47) becomes 

^ G " {bvvGy— X)' b, n G n J sin 0„ (A48) 

For symmetrically loaded wings, equation (A48) reduces to 

771+1 

A+ 


Gd { — , 1 Gm+l ( b^+l m +1 Gm+l — X+ -Bm+1 „ G„\- 

" l ~r 1 2 \ 2*2 2 n = l '2 ' / 


m+1 


'.'^ZGy( byyGy — '^,' By n G STO 0,1 
v=l \ n = \ / J 


(A49) 


Induced drag. — The induced drag is given by 

C d =A J* ^ G(jj)aiG])dri=A G(4>)<Xi(<i>) sin <f>d<t> 
or using the quadrature formula of equation (A15), 

<?z>‘ = m+I S Gya u sin 0„ (A47) 

For the straight wing, a t is the induced angle at the one- 
quarter-chord load vortex and the local-induced-dragdistribu- 
tion can be found. For the swept load vortex with the break 
at midspan it is not apparent how the correct induced angle 
at the load vortex can be found to determine local induced 
drag; however, the total induced drag of the wing can be 
found by considering the downwash in the wake of the wing. 

With the use of Munk’s stagger theorem (reference 16) 


where ^ are subscripts, and the prime on the summation 
sign indicates that the value for n=v is not summed. 

INTERPOLATION FUNCTION FOR SPANWISE LOADING 

The simultaneous solutions of the linear equations given 
by equation (A36) give the spanwise loading at wing 
semispan stations. The spanwise loading function is given 
by equation (A17) in terms of the known loading at — ^ 

points. An interpolation function can be determined to give 
values of loading between the known values and facilitate 
plotting of spanwise loading distribution. 

The spanwise loading given by equation (A 13) and the 
a H coefficients given by equation (A 16) can be arranged 
into the form 


GW 


2 

m + 1 


[(sin 0i sin 0+sin 20! sin 20 + . . . +sin m0! sin m0)G ! 1 + 


(sin 02 sin 0 + sin 202 sin 20+ . . . +sin w.0 2 sin m0)(?2+ 


(sin 0 m sin 0+sin 20 m sin 20+ . . . +sin m0 m sin m0) G m ] 


(A50) 


The values within the brackets can be tabulated for a given 
0 an d the loading at 0 will be the sum of the products of the 
tabulated constants and the known values of G n . 

For symmetric loading where only odd values of mi are 
needed equation (A50) becomes with m= 7 

(7(0)= i (0.383 sin 0+0.924 sin 30+0.924 sin 50+ 

0.383 sin 70) G^+i (0.707 sin 0+0.707 sin 30— 

0.707 sin 50-0.707 sin 70) <? 2 +-| (0.924 sin 0- 

0.383 sin 30—0.383 sin 50+0.924 sin 70) #3+ 

— (sin 0 — sin 30+sin 50 — sin 70) G t (A51) 

- Jctt 13 5 7 

Letting 0=0*=- g- for k= ^ and the factors of G„ of 

equation (A51) can be tabulated. 


The interpolation factor, e nlc , with m = 7 for symmetric 
loading is shown by the following table: 


V 

0. 981 

0.831 

0. 556 

0. 195 

\ k 





\ 

n \ 

\ 


3. 

2 

Jl 

2 

7 

2 

1 

0. 806 

0.490 

-0. 098 

0. 023 

2 

—.375 

.768 

.513 

-.075 

3. 

.278 

-.352 

.791 

.416 

4 

.128 

. 150 

-.225 

.641 


For symmetric loading 


4 

G[: ■ ^ ny,(+ n 


71 — 1 


(A52) 


Any type of loading coefficient other than G can be used 


equation (A52) , such as 


c ; c 
Glfiav 


in 
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For basic loading (C L = 0 ), equation (A46) gives 

£ 4 = -2(0.383 &+0.707 G 2 + 0.924 <? 3 ) (A53) 

With equation (A53), G t can be eliminated in equation (A51) 
and a shortened interpolation table for basic loading obtained. 

The interpolation factor, e„*, with m= 7 for symmetric 
basic loading is shown by the following table: 


V 

0. 081 

0.831 

0. 556 

0.195 

\ * 
n 

i 

3 

S 


_ 7 _ 

2 

1 

0. 963 

0. 375 

0. 075 

- 0. 468 

2 

-.195 

.556 

.831 

-.981 

3 

.513 

-.630 

1.207 

-.768 


For symmetric basic loading 

(A 54) 

71—1 

REFERENCES 

1. Weissinger, J.: The Lift Distribution of Swept-Baek Wings. 

NACA TM 1120, 1947. 

2. Van Dorn, Nicholas H., and DeYoung, John: A Comparison of 

Three Theoretical Methods of Calculating Span Load Distribu- 
tion on Swept Wings. NACA TN 1476, 1947. 

3. DeYoung, John: Theoretical Additional Span Loading Character- 

istics of Wings With Arbitrary Sweep, Aspect Ratio, and Taper 
Ratio. NACA TN 1491, 1947. 


4. Stevens, Victor I.: Theoretical Basic Span Loading Character- 

istics of Wings With Arbitrary Sweep, Aspect Ratio, and Taper 
Ratio. NACA TN 1772, 1948. 

5. Multhopp, H.: Die Anwendung der Tragfliigeltheorie auf Fragen 

der Flugmechanik. Bericbt S2 der Lilienthal — Ges. f. Luftf.- 
Forschg., Preisausschreiben, 1938-39, pp. 53-64. 

6. Margenau, H., and Murphy, G. M.: The Mathematics of Physics 

and Chemistry, D. Van Nostrand Co., Inc., New York, 1943. 

7. Anderson, Raymond F.: Determination of the Characteristics of 

Tapered Wings. NACA Rep. 572, 1936. 

8. Krienes, Klaus: The Elliptic Wing Based on the Potential Theory. 

NACA TM 971, 1941. 

9. Swanson, R. S., and Priddy, E. L.: Lifting-Suriace-Theory Values 

of the Damping in Roll and of the Parameter Used in Estimating 
Aileron Stick Forces. NACA ARR L5F23, 1945. 

10. Jones, R. T.: Properties of Low- Aspect-Ratio Pointed Wings at 

Speeds Below and Above the Speed of Sound. NACA Rep. 
835, 1946. 

11. Zimmerman, C. H.: Characteristics of Clark Y Airfoils of Small 

Aspect Ratios. NACA Rep. 431, 1932. 

12. Lange and Wacke: Test Report on Three and Six-Component 

Measurements on a Series of Tapered Wings of Small Aspect 
Ratio (Partial report: Triangular Wing). NACA TM 1176, 
1948. 

13. Shortal, J. A., and Maggin, B.: Effect of Sweepback and Aspect 

Ratio on Longitudinal Stability Characteristics of Wings at 
Low Speeds. NACA TN 1093, 1946. 

14. Glauert, H.: The Elements of Aerofoil and Airscrew Theory. 

The MacMillan Co., New York, 1943. 

15. Multhopp, H.: Die Berechnung der Auftriebsverteilung von 

Tragfliigeln. Luftf.-Forschg., Bd. 15, 1938. 

16. Munk, Max M.: The Minimum Induced Drag of Aerofoils. 

NACA Rep. 121, 1921. 





22 










24 














28 









30 











2 


,g a.o 

(m) »-4, n=*l,di=0.320. 
Figure 1. — Continued. 
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(n) v=i, n~ 2, d4=0.320. 
Figure 1— Continued. 
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Aspect-ratio parameter. 


THEORETICAL SYMMETRIC SPAN LOADING AT SUBSONIC SPEEDS FOR WINGS HAVING ARBITRARY PLAN FORM 


Note'- Basic loading characteristics of shaded wings are discussed in text. 
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A = 0 
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I — , — : 3 

A = 1.5 


-45 -30 O 30 45 60 75 

Compressible sweep parameter, An, degrees 

Figure 2. — Range of plan forms included in the study of the characteristics of straight-tapered wings. 
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•'I C 

Figure 3.—' Variation of the spanwise loading coefficient — — ^ — with the compressible sweep parameter Ap, degrees, for various values of the aspect-ratio parameter — and taper ratio X. 
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c u c 



(d) 1,-0.924. 


Figub* 3 —Concluded. 
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Spanwiae loading coefficient for unit .twist, " f > P er ^*3 






Spanwise loading coefficient for unit lift coefficient. 
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O .5 1.0 1.5 2.0 2.5 3.0 3.5 <0 4.5 5.0 

Aspect ratio, A 


Figure 10 —Comparison of lift-curve slopes as obtained by the simplified lifting-surface 
method of the subject report with those obtained by two lifting-surface methods and by a 
theory for wings of zero aspect ratio. 



Figure ll.— Comparison of lift-curve slopes obtained experimentally for triangular and 
rectangular wings of various aspect ratios with those given by simplified lifting-surface 
theory of subject report. 
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Theoretical aerodynamic-center location, a.c., fraction of M.A.C. 
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Experimental aerodynamic -center location , <L. c., fraction of M.A.C. 

Figure 13— Comparison of aerodynamic-center locations obtained experimentally for a random group of wing plan forms with those given by the simplified lifting-surface theory. 




Spanwiae 



Figure 15.— The effect of aspect ratio on the basic loading of wings having a taper ratio of 0.5 







Spanwise loading coefficient for unit twist. 
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Spanwise station, rj, fraction of semispan 
1 16 —The effectof taper ratio on the basic loading of wings having an aspect ratio of 3.6 


./ .2 .3 .4 .5 .6 .7 .8 .9 ,1.0 

Spanwise station , 77 , fraction of semispan 

Fiqurk 17— The effect of taper ratio on the basic loading of wings having an aspect ratio of 0.0 
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Figure 19.— -The effect of plan form on the angle of zero lift. 




Spanwise loading coefficient for unit lift coefficient. 
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(a) A = -45». 

(b) A=0°. 

Figure 20 —The effect of taper ratio on the additional loading of wings having an aspect ratio 
of 3.0 and various values of sweep. 



(d) A =60°. 

Figure 20. — Concluded. 











Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 


x 

X 


L 

Y-—*Z 

Roll 

rV 

. u 

V 


Y 

Y 

Pitching 

M 

Z *X 

Pitch.. 


V 

q 

Normal 

Z 

Z 

Yawing 

N 

X > Y 

Yaw 


w 

r 


AbsoLufce coefficients of moment 

c-*L 

°‘~qbS m ~qcS 

(rolling) (pitching) 


N 
qbS 
(yawing) 


C n = 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D Diameter 

p Geometric pitch 

p/D Pitch ratio 

V' Inflow velocity 

V, Slipstream velocity 

T Thrust, absolute coefficient C T — pn ijj* 

Q, Torque, absolute coefficient ^Q~ p ^]yi 


P 


C, 


V 

n 


Power, absolute coefficient C P — 


P 

P n 3 D 5 


Speed-power coefficient = 



Efficiency 

Revolutions per second, rps 


Effective helix angle = tan 1 



5. NUMERICAL RELATIONS 

1 hp=76.04 kg-m/s=550 ft-lb/sec 1 lb=0.4536 kg 

1 metric horsepower= 0.9863 hp 1 kg=2.2046 lb 

1 mph=0.4470 mps .•••••* 1 mi= 1,609.35 m=5,280 ft 

1 mps=2.2369 mph 1 m= 3.2808 ft 








